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Abstract 

Morphogen transport is a biological process, occurring in the tissue of living organisms, 
which is a determining step in cell differentiation. We present rigorous analysis of a simple 
model of this process, which is a system coupling parabolic PDE with ODE. We prove 
existence and uniqueness of solutions for both stationary and evolution problems. Moreover 
we show that the solution converges exponentially to the equilibrium in C 1 ' 01 x C°' a topology. 
We prove all results for arbitrary dimension of the domain. Our results improve significantly 
previously known results for the same model in the case of one dimensional domain. 

AMS classification 35B40, 35Q92 
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1 Introduction 

Morphogen transport (MT) is a biological process occurring in the bodies of living organisms. It 
is known that certain proteins (ligands) act as the morphogen - a conceptually defined substance 
which is responsible for the development of the shape, size and other properties of the cells. 
According to the 'French flag model' of Wolpert [T5], morphogen molecules spread from a localized 
source through the tissue of newly born individuals and after some time form stable gradients 
of concentrations. Receptors, located on the surface of the cells, detect those gradients and 
pass to the kernels the information about levels of morphogen concentration. Then according 
to these information, certain mechanisms begin synthesis of proteins which finally results in cell 
differentiation and specialization. Although the role of morphogen gradient in gene expression 
seems to be widely accepted, the exact kinetic mechanism of its formation is still not known, 
(see [S],[TU] and 0). 

Recently various models consisting of PDE-ODE systems were proposed to explain MT. Those 
models assume that movement of morphogen molecules occurs by different types of diffusion or 
by chemotaxis in the extracellular medium. Reactions with receptors (reversible binding, tran- 
scytosis) and various possibilities of degradation and internalization (of morphogens, receptors, 
morphogen-receptor complexes) are also being considered (see [IT], [8], 0, [IB]). 
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For the case of morphogen Decapentaplegic (Dpp) acting in the wing disc of the Drosophila 
Melanogaster individuals, several models have been proposed in In this paper we will be 
concerned with model [LNWJ.B (Model B [TT] p786). In this model it is assumed that movement 
of morphogen molecules occurs by passive diffusion while being affected by reactions of reversible 
binding with receptors and degradation of morphogen-receptor complexes. Morphogen is being 
delivered to the system by secretion from a source localized on one of the boundaries of the 
domain Q c K. n , which represents a fragment of the wing tissue. In mathematical terms the 
model is a system of two differential equations (PDE+ODE equipped with initial and boundary 
conditions), governing time evolution of the concentrations of free morphogen and morphogen- 
receptor complexes. 

In case of ID domains a detailed mathematical analysis of this model was made in [TJ] and 

m- 

In |14) the case Q = (0, oo), with a nonlinear dynamic boundary condition at x = and vanishing 
boundary condition at x — ¥ oo is considered. Well-posedness and L p {Tl) convergence of the 
solution to unique steady state were proved. 

In [7] the case Vl = (0,1), with nonhomogeneous, constant Neumann condition at x = and 
homogeneous Dirichlet condition at x = 1 is analyzed. Finding Lyapunov functional allowed to 
prove well-posedness and Z>2(^) exponential convergence to the unique equilibrium, with rate \ 
expressed explicitly by the parameters of the model. 

The goal of this paper is to examine [LNW] .B in the [7] setting for bounded domains of arbitrary 
dimension n. Although n g {1, 2, 3} is, from the biological point of view, the only relevant case, 
we do not put this restriction on n (methods that we use do not depend on the dimension). 
Using fixed point theorem and monotonicity of the nonlinearity we prove that our model has a 
unique nonnegative steady state. Using theory of analytic semigroups and comparison principle 
arguments we show existence of classical global solutions. We check that Lyapunov functional, 
obtained in [7J, also works for arbitrary n and thanks to appropriate semigroup estimates and 
bootstrap arguments we improve the topology of the convergence to the equilibrium from L 2 x L 2 
to C > a x C°' a without losing the exponential rate X- 

2 The model 

We consider the system of differential equations governing the space and time evolution of the 
concentration of free morphogen / and concentration of bounded receptors s in an annular shape 
domain Q C M™. We assume that receptors are distributed uniformly in the tissue so after 
normalizing the total concentration of receptors (free+bounded) is equal to I. The model governs 
the following biological processes 

• Passive diffusion of morphogens in the extracellular medium. 

• Secretion of morphogens from the source on a subset Tat of <?f2. 

• Binding of morphogens to receptors. 

• Unbinding of morphogens from receptors. 

• Degradation of bounded morphogens. 
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We equip the model with initial conditions Iq, so and boundary conditions on Yjj , r^v - two 
disjoint parts of <9f2. On Tjv we consider nonhomogeneous, time independent, nonnegative Neu- 
mann condition (flow of morphogen into the domain) while on r jj we put homogeneous Dirichlet 
condition (far from the source of morphogen their impact on the whole process is negligible). 
After normalization we end up with the following model 

[LNW].B 

d t l - DM = 5s - 1(1 - s) , (t, x) G (0, oo) x Q 

d tS = -(5 + e)s + l(l- s) ,(t,x) G (0,oo) x fi 

-DV n l = -v , (t, x) G (0, oo) X 

1 = , (t, x) G (0, oo) X Tp 

1(0) = l ,xeQ 

s(0) = so , x G O 

where we denote the derivative in the direction of the outer normal vector to Tn by V n . 



3 Results 

In the whole paper we assume that 
Al n G N, p > n > 1. 

A2 51 C K™ is a bounded domain (open, connected) with (C 1 ' 1 ) boundary which consists of two 
disjoint parts: d£l = U IV- 

A30<y£ Wp 1_1/p (rV). 

A4 i 0) so G < l (x), < s (a:) < 1, for x G O; Z (a;) = s (x) = 0, for x G Ud. 

Under the above assumptions we first analyze the stationary problem and prove the following 

Theorem 1. [LNWJ.B has unique nonnegative steady state [loo,Soo), where 
< ^oo € Wp(Q) is the unique solution to 

-DAl 00 = - x d °° ,xen (la) 

o + e + loo 

-DV n l 00 = -v ,x£T N (lb) 

Joo=0 ,x£T D . (lc) 

and Sqo = loo/(e + 5 + loo). 

The proof of existence is based on maximal regularity for uniformly elliptic operators in Sobolev 
spaces, compact embedding, comparison principle and Schaudcr fixed point theorem. Uniqueness 
follows from monotonicity of the nonlinear part in (|la[) . 



We next turn to the evolution problem and establish its well-posedness. 
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Theorem 2. [LNWJ.B has unique solution (l,s) such that 

I -loot C([0, oo); W^(Q)) n C\(0, oo); H^(fi)) n C((0, oo); VK p 3 (r!)) (2a) 
s6C 1 ([0,oo);^(n)). (2b) 

Moreover for (t,x) G [0, oo) x f2 

0<l(t,x), 0<s(t,x)<l. (2c) 

Local existence and uniqueness are obtained by putting system [LNWJ.B into the semigroup 
framework and using general theory for abstract parabolic semilinear problems. Comparison 
principle allows us to deduce that (|2"cf is satisfied from which we get that our solution is global. 

We finally study the stability of the steady state and show that it attracts all trajectories with 
the uniform exponential rate. 

Theorem 3. There exists a positive constant C depending on lo, so, v, 6, e, D, fl,p such that for 
every t > 

\\l(t) - loch, P + IK*) - soclkp < Ce-^\ (3a) 
\\l(t) - looh, P < Cmax{l/Vt, l}e-^ 2 )*, (3b) 

where 

. r_. DXi(6 + e) e-i fn . 

and Ai is defined in Lemma\^ 

By extending Lyapunov functional (derived in [J for one dimensional interval) to the case of 
arbitrary dimension we obtain estimates on the distance between solution and steady state in 
L2 x L2 topology. Using regularising properties of the heat semigroup we next bootstrap the 
topology of convergence to W% x Wp . 

Remark 

Using embedding W^(fl) x W£(Si) C C^ a (n) x C°' Q (ft) valid for p > n, < a < 1 - n/p we 
obtain topology of convergence as claimed in the introduction. 

4 Notation, semigroup estimates, Gronwall inequality 

For ijelwe denote x V y := max{a;, y}, x A y :— min{a;, y}, x + :— x V 0, x_ := (—x) V and 
extend this notion to real valued functions. If (V, >) is partially ordered vector space we denote 
its positive cone by V+ := {v € V: v > 0}. 

We make standard convention that C denotes positive constant which may depend on a subset 
of {/q, sq, v, S, e, D, fl,p} and may change its value from line to line. 
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For 1 < q < oo, a G {1,2,3} we introduce the spaces W^ Ba (fi): 

{uGW^fi) :u|r D -0} 
{w G W g 2 (Q) : u|r D = 0, V„u|r N = 0} 

G W g 3 (Q) : u | Fd = 0, V n u\ TN = 0, A u |r D = 0}, 

with standard Sobolev norms ||.||a,g- 

We next recall some properties of the heat semigroup generated by laplacian with appropriate 
boundary conditions. 

Lemma 1. For 1 < q < oo the Laplace operator A: L q (Q) D W^ B2 (^) — > L q (Q) generates an 
analytic, strongly continuous semigroup e tA . For a, ft G {0,1,2,3}, a. < ft, 1 < qi < q2 < oo 
and t > we have 

||e tA «lk 9 < C(t A l)f a -«/ a e- Alt ||«|| 0l , < Ct^-^ 2 \\u\\ a , q , u G W£ fl<> (4a) 

||e tA u|| g2 < C(t A l)-" /2(1/ « 1 - 1/ « 2) e- Alt ||u|| gi < Ct-^ 1 /' 1 -^Hullft, u G L ft (4b) 

where X% > is the first eigenvalue of —A and C depends only on q, qi, 02, ^- 

Proof. Noticing that — Ai = sup Re(a(A)) we get from [TJ following estimates 

lie tA u|| 9 <Moe- Al *||M|| g 
||t(A + AiJ)e* A u|| ? < Mie -Al *||u||,, 

We have 

\\e tA u\\ 2 , q < C\\Ae tA u\\ q < C\\(A + XiI)e tA u\\ q + C\ x \\e tA u\\ q 
< C{Mx/t + MoAi)e -Al *||u||g < C(t A l)- 1 e _Al *||tt|jg 

From [1] we have that 

[L g ,Wl B2 ] a/2 = W« Ba , a G {0,1,2,3}, 
[L qi ,Wl tB2 ] e C L q2 , 9 > n/2(l/ qi - l/q 2 ), 

where for 9 G [0, 1] [.,.]g denotes complex interpolation functor, which is extended for 9 > 1 as 
described in [Tj. From this estimates (|4al) and (I4bp follows. □ 

We next recall the singular Gronwall inequality 

Lemma 2. Assume that f G C([0,T);]R + ) satisfies for every t G [0, T) following inequality 

f(t)<a + b [ (t-s)- a f(s)ds, 
Jo 

where a,b are nonnegative constants and a G [0,1). Then there exists positive constant C — 
C{h, a) such that for t G [0, T) 

u{t) < aCe bct . 

Moreover C(b,Q) = 1. 



w£ B i(n) = 
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Proof. For proof (under more general assumptions) see Lemma 7.1.1 in [BJ. 



□ 



5 Proof of Theorem [T] 

For x > let f(x) = s+ e c+x ■ Consider the operator T : L p (Q) + — > L p (Q), denned by T(u) = w 
where w € W p (f2) is the unique solution of 

-DAw + f(u)w = ,xeQ (5a) 

—D 1 \/ n w = —v ,x^Tn (5b) 

w = ,x£T D (5c) 

We will show that T has bounded range in L p (Q) + , is compact and continuous (this via the 
Schauder theorem will imply existence of a solution of ([T]) in Wp(£l)). Using the fact that 
< f(x) < we get from maximal regularity of uniformly elliptic differential operators in 
Sobolev spaces (see [4] for instance) the following estimate 

IMIwj(n) < C M w ^ 1/p (r N ) 

which gives boundedness of the range of T in W p (f2) and therefore in L p (fl). Compactness of T 
follows from the compact imbedding W p (Vl) CC L p (Q,). To show that w > we multiply (|5"aj) 
by u>_ and integrate by parts (notice that for p > n w £ W p (SY) C (fi) hence u>_ € WaO^O) 
to obtain 

-£> / |Vw_| 2 - / vw-- I f(u)w 2 _ = 0. 
Since w = on therefore u> > in Q. 

Assume that u„ — ► u in L p (Q). Let u> = T(u),w n — T(u n ), then 

-DA(w n -w) + f(u n )(w n -w)+ w(f(u n ) - f(u)) = , x € 

-DV„(w„-w)) = o ^erw 
w n — 1« = ,ierj) 

therefore 

||w„ - HU P (n) < C|K/(u n ) - f(u))\\ Lp{n) < C\\w\\ Loa{n) \\f'\\ Lx{0tOo) \\u n - u\\ Lp{n) 

which proves that T is continuous. Using Schauder fixed point theorem we obtain existence of 
loo € W p (£l) which solves §3). 

To prove uniqueness, assume that lloJlo are solutions of (P). Subtracting equations ([Taj) for 
^ooi^oo > multiplying by ^ — integrating by parts and using the monotonicity of function 
K + 3i-> xf(x) we get 

-d f |v(/L-4)l 2 = / (/(&)&- -&)> o, 

which by (ITcl) implies = 
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6 Proof of Theorem [2] 



To deal with nonhomogeneous boundary condition on Tjv we subtract from (I, s) the stationary 
state (Zoo, Soc). Setting (zx, z 2 ) = (I — loo, s — Sqo) we arrive at 



- DAz! 


= Sz 2 - 


2l(l - 2 2 ) + SooZl + ZooZ 2 


,(t,x) e (o,oo) 


x Q 


(6a) 


d t z 2 


= -(54 


- e)z 2 + 2l(l - Z 2 ) ~ Soo2l ~ Zoo22 


,{t,x) G (0,oo) 


x n 


(6b) 


-DV n Zl 


= 




,(t,x) G (0,oo) 


x r N 


(6c) 


Zl 


= 




,{t,x) G (0,oo) 


x r D 


(6d) 


zi(0) 


= 210 = 


^0 — ^oo 


,x eft 




(6e) 


Z2(0) 


= 220 = 


SO - Sqo 


,x e 




(6f) 



We interpret system ([5]) as a differential equation in a Banach space specified below 

z-Az = H(z) ,*e(0,oo) (7a) 

z(0) = z = (210,220) (7b) 

where z = (zx, z 2 ), Az = (DA Zl , 0), H = (H 1 ,H 2 ), 

H l {z) = Sz 2 - zi(l - z 2 ) + SooZi + Zoo22 (8a) 
H 2 (z) = -(5 + e)z 2 + 21(1 - z 2 ) - S0021 - Zoo2 2 . (8b) 

In the following lemma we prove local existence for (J7]). 

Lemma 3. For a G {0,1,2,3} denote Z a _ p = W^ B ^ x W^ Bl . For every z G Z\ tP the Cauchy 
problem ([7]) possess a unique maximal local solution 

z G C([0, T max ); Z hp ) n ^((0, T max ); Z liP ) n C((0 

, -^iriax)) ^3,p)* 

which satisfies for t G [0, T max ) the following Duhamel formula: 

zi(t) - e* DA z 10 + / e( t -^ DA i/ 1 (z(s))ds (9a) 
io 

22(f) = 2 20 + / H 2 (z(s))ds. (9b) 
./o 

Moreover «/T max < 00 £/ien limsup i _ s . T - ||z(f)||i,p = 00. 

Proof. The operator A : Z p Z) Z 2tP — > Z p is a generator of an analytic strongly continuous 
semigroup e tA — e tDA x Id (as a product of two generators). Moreover since Z\ jP is a Banach 
algebra (p > n) we observe that i7 : Z\ tP — > is locally Lipschitz on bounded sets. The claim 
follows from Theorem 7.2.1 in j3J. □ 

We next turn to the proof of (|2"c|) . 
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To prove that for t G [0, T max ) l(t), s(t) > we consider the system 

d t l' -DAI' = Ss' + -l' + (l-s' + ) , (t, x) G (0, oo) x fl (10a) 

d t s' = -(5+e)s' + + l' + {l-s' + ) ,(t,x) G (0,oo) xfi (10b) 

-DVJ' = -i> , (t, x) G (0, oo) X T N (10c) 

/' = G (0,oo) x Tz, (lOd) 

Z'(0) = Z ,xefi (10e) 

s'(0)=s ,2; eft (lOf) 

As before one can show that (|10[) possess unique classical local solution (/', s'). After multiplying 
(jlOaj) by Z_ and integrating by parts we obtain 

-~ [ \l'_\ 2 dx-D [ \Vl'_\ 2 dx-D [ l'_vdS = 8 [ s' + l'_dx>0. 
2 dt J n Jq J Tn J n 

Similarly multiplying (jlObp by s_ yields 

1 d 

"Ydt 



s'_\ 2 dx = I l' + s'_dx > 0. 

Q 



Therefore for t G [0, T max ) 

Iin*)-lla+ll*'(*)-ll2<l|i'(0)-||i + ||*'(0)-||i = 

and consequently l'(t) > 0, s'(i) > 0. We observe now that (V, s') is a solution of [LNWJ.B and 
using uniqueness we finally get that l(t) = l'{t) > 0, s(t) = s'(t) > for t € [0, T max ). 

To show that s(t,x) < 1 for (t,x) G [0,T max ) x we get from Lemma 02 that for every fixed 
x € fl the function s = 1 — s = 1 — Z2 — Soo G C 1 ([0; T max ),R) satisfies for t > the following 
ODE 

s + (5 + e + l)s = S + e. 

Therefore 

S(t) = e-( 5 + £ )*-/o 'W<Jr (1 _ Sq) + (J + e ) /** e -(i+«)(*-*')-; «- t ' Kr)dT dt , > Q 



We finally show that T max = oo. Reasoning by contradiction assume that T max < oo. Using 
uniform boundedness of s (and therefore of z%) we obtain for t G (0,T max ): 

II^W*))!!, < c(i + \\ Zl (t)\\ p ) < c(i + lkiWlli,p)- (ii) 

Using (|9a). (|4a) .([TTj) we obtain 



lki(*)lli, P < l|e* DA ^io||i. P + / WeWWMTMxjdT 

Jo 

<C\\zioh, p + C f (t-Ty^WH^zmipdr 



< C\\z w \\ liP + C / (t-T)- 1 / 2 {l + \\z 1 {T)\\ 1 ^)dT 

Jo 

< C(||2io||l,p + 1) + C [ (t - T)-^ 2 \\zi{T)\\ lip dT 
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Using Lemma [5] we get that < C and therefore 

||^ 2 (^))||i >P <C(l + || 22 (t)|| 1>p ). (12) 

Using (J9bj) and JT2J we obtain 



\\Z2{t)\\l >p < \\Z20\\l, P + / \\H*(z(T))\\ llP dT<\\z 20 \\l, P + C (1 + ||z 2 (r)|| l, p )dT 

Jo Jo 

<C(\\Z 20 \\ 1 . P + 1)+C f \\ Z2 (T)\\ hp dT. 

Jo 

Another application of Lemma [5] gives desired contradiction from which we deduce that T n 
oo. 



7 Proof of theorem [3] 

The proof of Theorem [3] is based on L 2 estimates obtained for n ~ 1 in [7] and bootstrap method 
to improve convergence from A^-topology to Xi + i-topology, where Xi + \ c Xi are appropriately 
chosen Banach spaces. We use (as long as the regularity of our solution permits) the following 
two step 

Bootstrap scheme 

1. + \\z2{t)\\ Xi < Ce-WW gives \\ Zl (t)\\ x . i+1 < Ce~^. 

2. |M*)lk +1 < Ce-WW gives \\z 2 (t)\\ Xi+1 < Ce-^K 

Step 1. is a consequence of Duhamel formula (|9"aj) and semigroup estimates (j3]). 

Step 2. follows from the fact that we can solve equation (|6b[) explicitly for z 2 in terms of Z\. 

7.1 L 2 estimate 

We first show that, as in the one dimensional case [LNWJ.B has a Lyapunov functional from 
which exponential convergence to the equlibrium (^oo,Soo) follows. 

Lemma 4. For x € [0, 1), u, v £ W p x B i(fi), < v < I, define 
S 7 (ic) = -ln(l - a;) 

r i) — "I 

dx 



1 - Sr 



Ao(«)= / (l-s 00 )(/ 00 + 5 + 2e) S 7 ( v )~E / ( Soo )- W 
Jn 

A(u,v) = t;\\ u - l °o\\l + A o(«) 



T>k{U,V) = D \\\{u- loo) 2 + / 1 u^- 



T/ien /or t > 



A(J(t), «(*)) + / Pa(Z(t),s(t))^ = A(Zo,so) 



X A(i(i),*(t))<I> A (l(i),s(t)) 
(5 + e)||s(t) - Soo ||l < 2A (s(t)) 
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and 

WKt) - UI2 + (<* + ~ Soolls < 2AG , s )e-*\ (13) 

where \ satisfies (|3cl) . 



Proof. Proof can be obtained exactly as in [7] (part of Theorem 8 and Proposition 9 pp 1740- 
1744). For the case n = l,p G (1, 2), to justify integration by parts and Poincare inequality, we 
observe that for t > : l(t) G W%{Sl) C W£(Q). □ 

7.2 L p estimate 

In this subsection we will prove that for t > 

ll*i(*)l| P +ll«a(*)llp<C'e-(x/ 2 ) t ) (14) 
the parameter p being defined in Al. 

Notice that if p e (1,2] (which can only happen if n = 1), the inequality (|14p follows from 

JH. 



Otherwise we have p > (2 V n). We choose an increasing sequence (pi)i^i such that 

Pi = 2,p m = p 
n/2(l/ Pi - 1/pi+i) < 1 

(notice that for n G {1, 2, 3, 4} one can take m = 2). Inductively we will prove that 

ll«i(*)ll« + \\z2(t)\\ Pi < Ce~W)\ 1 < i < m. (15) 
For i = 1 (fT5"j) follows from (fTl?| . Assume that (fTS")) is true for some 1 < i < m — 1. Then 

\\RH*(t))\\pi < INUI|l-«2 + floo||oo + INUH« + e + ioo||oo < Ce-^ 2 )*. (16) 
Using ([5a]). (|4b|) . and x/2 < DX 1 we obtain 

IM«)|| P<+1 < \\e tDA z 10 \\ Pi+1 + fv D *H\z{t-s))\\ Vi+1 ds 

Jo 

<Ce- DXlt + C f {Ds Aiy n ^ 1 ^'~ 1 / p '^e^ DXlS \\H 1 (z(t - s))\\ Pz ds 
Jo 

< Ce- DXlt + C f (Ds A iyn/WPi-i/Pi + i) e -DX 1Se -( x /W-s) ds 

Jo 

< Ce- DXlt + Ce-W>* [\ds A i)-"/2(i/ft-i/p (+1 ) e -(i3A 1 - x /2)^ s 







To show that for t > ||z2(i)|| Pi+1 < Ce (x/ 2 )*, we obtain from Theorem [2] that for each fixed 
x G f2 the function z 2 G C x ([0, 00); R) satisfies the ODE 

z 2 + (5 + e + loo + zi)z 2 = (1 - 800)2:1, 
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hence 



Z 2 (t) = A(t) Z2Q + (1 - Soc ) f A(T) Zl {t-T)dT, (17) 

Jo 



where 

A(t) = exp ( - J (6 + € + 1^ + zi(r))dr) . 
From Zoo + zi = I > we get U^i)^ < e -(<S+e)t, Using x/2 < <5 + e we obtain 



(18) 



H^WHw+x < ll^(*)l|oo||«20|| W+1 + ill - SooHoo / ||A(T)|U|zi(t - T)|| pi+1 dT 

< Ce -( 5 + e )* + Ce-^ 4 f e-( 5+e -x/2)r dT < Ce -(xmt 



ii 



thus finishing the proof of (fI5]l. whence that of (JT? 



In the next two sections we use the smoothing properties of e tA to extend convergence to the 
first and second derivatives. 

7.3 Wp estimate 

Using ([5a]l. gaj), flT3J| and x/2 < DAi we obtain 

rt 

JBA. , / llsDA rrl/ 



lki(t)||i )P < II enroll r, P + / \\e s ^H'(z(t-s))\\ hp ds 

Jo 

< Ce- DX ^ + C [ (Ds A lr^e-^Hff 1 ^ - «))|| p dfl 

Jo 

< Ce~ mit + C ['(Ds A ij-i/ae-Ai^e-^Ct-.)^ 



< C-g-DAit + Ce -( X /2)t / ( DsA1 )-V* e -(D^-Xms ds 







Using the above estimate for z\ we obtain that A(t) given by (fT5)l satisfies 
II^WIIp < CptoUoo < C e -<*+ e >* 

||V4(t)||p = ||-A(t) / (V/oo + V2l(T))dT|| p < IIAWHoo / (HVIoollp + ||V2l(T)|| p )dT 

Jo Jo 

< Ce"^* Al + e^^dr < ffte"^'. 
Jo 
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Thus using (jTTJ) we have 

IM*)||l,p < m(*)||l I p||220|K,p + - SooHi.p f ||A(r)||i iP ||^i(t - T)\\l, p dT 

Jo 

< C(t + l)e-( 5+£ )* + C [\t + l^+^e-^^dT 

Jo 

< C(t + l)e-( 5+e )* + Ce-^ 1 [\t + l)e- {&+ '-^ T dT 

Jo 

< Ce-^ 1 
which finishes the proof of (l3al) . 

7.4 Wp estimate for z\ 

Using ([5a]). (gaj), ([21]) and x/2 < L>Ai we obtain 

lki(*)||a,p< \\e tDA z lQ \\ 2 , p + I ||e TDA H 1 (2(t-T))|| a , p dT 

Jo 

< C7(Dt A l)-l/2 e -BAlt + C A ir l/2 e -A 1 I3r e -( X /2)( f -r) ds 

JO 

< C(i A l)-V2 e -SA lt + Ce -(x/2)t f*( T A 1) -l/2 e -(DA 1 - X /2)r dr 

Jo 

< Cmax{l/Vi,l}e" (x/2) *, 
which finishes the proof of (13b|) . 
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